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Abstract. We investigate the motion of free relativistic particles relative to the ambient medium around a gravitationally 
collapsed system. If the relative speed exceeds a critical value given by c/ V2, the gravitational tidal effects exhibit novel 
features that are contrary to Newtonian expectations. In particular, ultrarelativistic jet clumps moving freely outward along 
the rotation axis strongly decelerate with respect to the ambient medium, while ultrarelativistic particles strongly accelerate in 
directions normal to the jet axis. The implications of these direct consequences of general relativity for jets in microquasars and 
the origin of the high-energy cosmic rays are briefly mentioned. 
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The purpose of this Letter is to discuss some of the observ- 
able consequences of the motion of relativistic particles in the 
field of gravitationally collapsed configurations. In particular, 
we point out novel features of general relativity regarding the 
tidal influence of a neutron star or a black hole on the motion 
of nearby free particles moving faster than the critical speed 
c/ V2 ss 0.7c relative to the ambient medium surrounding the 
central mass. Such ultrarelativistic particles originating near 
the poles of the collapsed system suffer significant tidal de- 
celeration within a cone of angle a; 55° around the axis of 
rotation. Here 6 represents half the angle at the vertex of the 
cone. Outside the cone, however, an ultrarelativistic particle ex- 
periences tidal acceleration. The situation is depicted schemat- 
ically in Fig.[2for the north pole of a Kerr black hole. Similar 
outflows are expected near the south pole. 

The solutions of the geodesic equation as well as the 
Lorentz force equation in black hole fields in standard 
Schwarzschild-like systems of coordinates have been exten- 
sively studied and many useful results have thereby been 
obtained (see, e.g., Chandrasekhar J983l. For the study of 
certain phenomena, however, it is more useful to study relative 
motion, i.e. motion of one particle with respect to another. 
This is in keeping with the spirit of the theory of relativity. 
It also corresponds to the physical interpretation of some 
astronomical situations. Consider, for instance, the Chandra 
X-ray images of the Crab Nebula: to reveal the accelerating 
rings in the equatorial plane, different images have to be com- 
pared with each other keeping certain central features "fixed" 
(http://chandra.harvard.edu/photo/2002/0052/index.htmli. To 
study relative motion in an invariant setting corresponding to 
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Fig. 1. Schematic diagram of the critical velocity cone. 



actual observations, it is useful to establish a Fermi normal 
coordinate system along a reference geodesic O. The motion 
of a nearby free particle with respect to O is then given by 
the geodesic equation in the quasi-inertial Fermi coordinate 
system (Synge 1960). 



Let T be the proper time and A 



(0) 



dx^/dT be the four- 



velocity vector of O. A triad of ideal orthonormal gyro direc- 



tions /l'',,,, / = 1,2,3, can be parallel propagated along <9's 
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worldline such that A'^.Jt) is a local orthonormal tetrad frame. 
Based on these local axes, the Fermi coordinates (T, X) simply 
provide a geodesic normal coordinate system along the world- 
line of O. The reference worldline is then the temporal axis 
(T - T, X - 0) and the metric tensor in the Fermi system is 
given by 

^00 = -l-%mj(T)rx^ + ---, 

goi = --%)jdT)XJx'' + ---, 



?U = Sij--%kjc{T)xV + ' 



where 



RaBysiT) - Ruypa-A'' A" ,g.A^ A"^ 



pypo-'^ (a/^ (fiY^ ^yy^ {6) 



(1) 



(2) 



is the projection of the Riemann tensor on the local frame of 
O. Along (3's worldline, the metric reduces to the Minkowski 
metric according to O. The deviations from an inertial system 
of coordinates are quadratic and higher order in distance away 
from the reference worldline. The Fermi coordinates are ad- 
missible only in a cylindrical region of radius 1^(7) around the 
reference worldline, where "R is a certain minimum radius of 
curvature determined via the Riemann tensor and its covariant 
derivatives along the worldline of O. The equations of general 
relativity may be expressed in any admissible coordinate sys- 
tem. We are interested in the geodesic equation of motion in 
the quasi-inertial Fermi coordinate patch. 



evolution of the deviation between two neighboring geodesic 
worldlines with negligible relative motion. Equation (|5} has 
been employed extensively in the treatment of tidal effects 
in general relativity (e.g. Stewart & Walker 1973, Mashhoon 
[T975l[T977l Carter & Luminet [T985l Ivanov & Novikov|200D- 
Detailed calculations have revealed that the dominant tidal ef- 
fects of massive bodies are rather similar to their Newtonian 
counterparts except possibly very close to black holes. 

Equation (|4} follows from the time like character of particle 
motion; that is, the requirement that g^yU^U^ - -I leads to (0) 
when we use the terms in the expansion of the metric tensor 
given explicitly in Q. The modified Lorentz factor is positive 
and approaches infinity if the timelike geodesic approaches a 
null geodesic. It follows from (0J that V^ < 1 along the refer- 
ence geodesic; elsewhere in the Fermi frame, however, V^ is 
simply constrained by the requirement that the right-hand side 
of (0} be positive. 

The quasi-inertial Fermi coordinate system has been widely 
used in the theory of general relativity; for example, it is em- 



ployed by Ehlers & Rindler (I1997> to study the local bending 
of light in a gravitational field. 

3. Tidal acceleration in Kerr spacetime 

In this Letter we point out the major consequences of (|3j and 
(|4} for ultrarelativistic motion in the field of a rotating source 
represented by the exterior Kerr spacetime. The Kerr metric can 
be expressed in Boyer-Lindquist coordinates {t, r, §, (p) as 



2. Generalized Jacobi equation 

The geodesic equation for a free test particle V in the Fermi 
system can be reduced to the general tidal equation (Mashhoon 
119751 11977> ; however, if we limit our considerations to the 
terms given explicitly in equation Q, we get the generalized 
Jacobi equation (Chicone & Mashhoon^2002|) 

-^ + ^RoiOjXJ + 2 %kjoV''X^ 



+ t(3 RokjQ 



V'V 



%kjcV'v' + %>,jerv'v')XJ = 0, (3) 



where V - dX/dT is such that the four-velocity of "P is given 
by V = r(l, V) and F is the modified Lorentz factor given by 



^- 



V' 



'RoiOjX'XJ + -'RojikX^V'X" 



+ -''RikjeV'x'VJx' 



(4) 



at the level of approximation under consideration here. For 
|V| <K 1, we may neglect all terms proportional to the rela- 
tive velocity in the equations of motion and Q reduces to the 
famiUar Jacobi equation (SyngeQ^^ 



^.^,,x^ = o. 



(5) 



Here Kjj = '^RoiOj is the symmetric tidal matrix that reduces in 
the Newtonian limit to d^O/dX'dX^, where <1> is the Newtonian 
gravitational potential. The Jacobi equation expresses the linear 



+ 2GM-{dt - a sirv^ &d(f>f. 



(6) 



where S - r^+a^ cos-^ -9- and A - p--2GMr+a^. Here Mis the 
mass, J - Ma is the angular momentum of the source and we 
use units such that c - \. We take (3 to be a test particle moving 
along the axis of rotation {§ - 0) on an escape trajectory such 
that it reaches infinity with a speed that is small compared to 
unity; in fact, we set this speed equal to zero, since our results 
turn out to be insensitive to it. The worldline of the reference 
particle is thus characterized by 



dt 



p- +a^ 



dr 



dT r^ - 2GMr -H a^ ' dr 



2GMr 



(7) 



We assume that at f = t = 0, the reference particle starts from 
r = ro > V3 a along the rotation axis and eventually escapes 
to infinity in the absence of other interactions. Given this refer- 
ence trajectory, it is then necessary to choose the spatial axes of 
the tetrad frame A'', ,. Taking A'',.,, to be along the direction of 

(a) o (3) o 

motion, axial symmetry about this direction implies that there 
is rotational degeneracy in the choice of A^ ^. and A''.^.. Once 
a triad is chosen, however, it is parallel propagated along the 
reference worldline. The projection of the Riemann tensor on 
this tetrad, as in (|5}, is given by a symmetric 6x6 matrix with 
indices that range over the set {01,02,03,23,31, 12): 



E B 
B -E 



(8) 
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where E and B are symmetric and traceless 3x3 matrices 
that represent the electric and magnetic parts of the curvature 
tensor. It turns out (Mashhoon, McClune & Ouevedo il999i ) 
that E and B are both diagonal and parallel; Ejk — B/q - 
diag(-5, -5,1), where 



-2GM 



r(r^ - 3a2) 



(r2 + a2) 



i\i 



q^l GMa 



3r2 



(r^ + fl2) 



2\i- 



(9) 



Integrating the equation of motion of the reference 
geodesic 0, we find r{T) and its substitution in equa- 
tion (|9jl gives the electric and magnetic curvatures k{T) < 
and q(T) > 0, respectively. The class of such reference 
particles provides an ambient medium around the source. In 
effect, this medium provides a local dynamic connection with 
the rest frame of the source. 

The generalized Jacobi equation Q in the case under con- 
sideration takes the form 



X--kX 

2 



1 



lie- + -(2Y^ 
3 



■z') 



+ -kX(5YY ■ 



7ZZ) 



+q[XYZX - ZY{\ + X^) - 2YZ] = 0, 



(10) 



Y--kY 

2 



1 



■ 2Y^ + -{2X^ ■ 
3^ 



Z^) 



1 



+ -kY(5XX - IZZ) 



-q[XYZY -ZX(l + Y^) - 2XZ] = 0, 



Z + kZ 



1 



2Z^ + -iX^ + Y^) 
3 ' 



H- -kZ(XX + YY) 



-q(XY -XY){\ -Z0 = 0, 



(11) 



(12) 



where X - dX/dT, etc. Let us first concentrate on the motion 
of !P along the rotation axis; that is, X{T) - Y{T) - for all 
r > 0. In this case, equations J10l l-J12> reduce to 






2Z})Z ■■ 



0. 



(13) 



The character of the solutions of this equation depends dras- 
tically on whether the relative initial velocity along the Z- 
direction has magnitude below or above the critical speed 
II yl. For initially infrarelativistic motion (i.e. Z} < 1/2), 
equation \\3\ implies that the particle accelerates until it 
asymptotically approaches the critical speed. For initially ul- 
trarelativistic motion (i.e. Z} > 1/2), however, the particle 
strongly decelerates regardless of its initial Lorentzfactor Fo > 
yl2 and asymptotically approaches the critical speed 11 \2. 
This deceleration phenomenon is consistent with the observa- 
tions of Galactic "superluminal"jets (Fender i2004i i. 

Let us next consider motion along the X-axis; that is, 
Y{T) = Z{T) = for all T > 0. Then equations ([Tol-(tT2i 
reduce to 



fx_ 

dT^ 



\kil 



2X^)X : 



0, 



(14) 



which turns out to be a generic equation for motion perpendic- 
ular to the symmetry axis as a consequence of the rotational 
symmetry of equations J10 M12> about the Z-axis. For initially 
infrarelativistic motion, it follows from \\A\ that the particle 
slowly decelerates, while for initially ultrarelativistic motion, 
the particle strongly accelerates as depicted in Fig.|2 The latter 
situation is remarkable, as it can lead to F — > 00 in the absence 




Fig. 2. Plot of the Lorentz factor F = {\-X^-kX^ I2)-^I^ versus 
T/(GM) based on integration of equation (I14> with initial data 
X = at r = with X(0) = (F^ - ly/^/To corresponding 
respectively to F,, = 5, 10, 15, 20, 25, 40, 70 and 100. In this 
plot a/(GM) - 1 and ro/(GM) - 15. The graph illustrates 
acceleration of the particle such that F essentially approaches 
infinity at T/(GM) * 1 1, 9, 7, 5, 3 and 2.6 for Fq = 15, 20, 25, 
40, 70 and 100, respectively. 



of other interactions. It is clear from Fig.|2lthat in some cases 
the Fermi coordinate system could break down before F — > 00. 
The breakdown in our dynamics for F — > 00 is logically distinct 
from the kinematic breakdown of the Fermi coordinate system 
beyond |X| ~ ^. These tentative conclusions follow from our 
simple scenario that takes into account only the quadratic terms 
in equation Q; in principle, the whole infinite series must be 
considered. That is, the geodesic equation of motion of T in- 
volves an infinite series of tidal terms; the breakdown occurs 
when instead of the infinite set of tidal terms only the first- 
order terms are taken into account. Our analytical and numer- 
ical studies of the general structure of this infinite series have 
led us to the working hypothesis that for ultrarelativistic parti- 
cles, finite but significant tidal acceleration takes place normal 
to the jet direction (Chicone & Mashhoon l2005> . 

The breakdown of our simple theoretical scheme as F — > 
00 implies that the influence of such a particle on the back- 
ground spacetime cannot be neglected. A more complete treat- 
ment should take this influence into account, thereby moder- 
ating the singularity that we have encountered. Moreover, the 
tidal energy exchange resulting in the deceleration/acceleration 
phenomena must be augmented to take due account of grav- 
itational radiation energy emitted by the particles. This en- 
ergy, as well as the corresponding electromagnetic radiation 
energy emitted by the electrically charged particles, is in the 
form of long-wavelength radiation of frequency proportional to 
(2GM)-' = 10^(Mq/M)Hz. The tidal acceleration of charged 
particles could be strongly affected by the electromagnetic field 
configuration around the source. 

It is interesting to consider the exchange of energy between 
the particle !P and the gravitational field. Einstein's principle 
of equivalence prevents a local (i.e. pointwise) description of 
gravitational energy. From the standpoint of observers at rest 
with the ambient medium, however, an initially ultrarelativistic 
particle would lose energy to the gravitational field along the 
jet direction, but would gain tidal energy propagating outward 
in the equatorial plane; thus, the gain in gravitational energy by 
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the accelerating particles could be essentially balanced with the 
loss of gravitational energy by the particles decelerating along 
the jets. 

For an invariant characterization of these tidal phenomena, 
consider the class of fundamental observers S that are at rest 
with four-velocity m^ - (W, 0, 0, 0) with respect to the Fermi 
coordinate system. Here W - (-goo)"'^^. Let f := —u^Ufj be 
the Lorentz factor of V measured by the static observers S, 
which are in general accelerated. It turns out that 



f = (W-' - WgQiV')Y. 



(15) 



We find that f behaves essentially as F; e.g., in Fig.|2 the plot 
of r would be indistinguishable from that of F. 

Ultrarelativistic particles are expected to be produced near 
a highly magnetized rapidly rotating neutron star or as a conse- 
quence of the complicated accretion phenomena in the vicinity 
of an active black hole (Punsly|2001 Guthmann et al. 2002i. 
For our purposes, we imagine an abundance of such particles 
near the poles of the collapsed system. Extensive numerical 
studies of equations J10t -J12> indicate that the deceleration 
along the rotation axis extends to a cone with an opening half- 
angle 6 given by tan0 a V2; moreover, motion outside the 
cone corresponds to tidal acceleration away from the collapsed 
system. 

Tidal acceleration of ultrarelativistic particles is maximum 
in the equatorial plane and decreases away from it, turning 
to deceleration at an inclination angle of about 35°. Our 
results regarding tidal acceleration appear to be consistent 
with recent Chandra X-ray studies of four neutron stars: 
Crab Pulsar, Vela Pulsar, PSR B 1509-58 and SNR G54. 1+0.3 
( http://chandra.harvard.edu/photo/chronological.html). In 

each case, there is considerable activity near the equa- 
torial plane of the central source. Moreover, a detailed 
analysis of the Chandra X-ray images of the Crab 
Nebula suggests that the equatorial activity is somewhat 
more energetic than the activity along the jets (Mori et 
al. 2004), which is consistent with the tidal accelera- 
tion/deceleration of ultrarelativistic particles studied in this 
work. The same appears to be the case for the Vela Pulsar 
( http ://chandra.harvard. edu/photo/2000/vela/index.html i . 

4. Discussion 

It is important to emphasize that the deceleration/acceleration 
phenomena are generic to collapsed configurations; the rota- 
tion of the central source has a negligible influence on these 
results. In particular, the results presented in Fig. |2 are insen- 
sitive to the assumption that the two particles initially coin- 
cide, or that the central source rotates rapidly. Numerical ex- 
periments for initially nearby particles as well as for a in the 
range [0, GM} demonstrate that the results presented here are 
robust. What matters is that the central source is a neutron 
star or a black hole; the largest effects are expected to occur 
in the latter case. In either case, tidal deceleration/acceleration 
can be significant near the source, where the ultrarelativistic 
particles are launched; however, a complete treatment should 
take due account of the complicated plasma processes as well 
(Punslv l200n Guthmann et al. l2002 i. 



The deceleration/acceleration phenomena occur relative to 
the ambient medium. That is, the geodesic orbits of the exterior 
Kerr spacetime under consideration in this work would not ap- 
pear to have any extraordinary features when considered from 
the standpoint of the static inertial observers at spatial infinity. 
The situation is different, however, in the Fermi system. 

The ultrarelativistic particles that are created near a 
gravitationally collapsed system or in the accretion pro- 
cess around the system are tidally decelerated in a cone 
around the rotation axis of the collapsed system and 
appear as confined relativistic outflows from the sys- 
tem. These correspond to jets from some neutron stars 
(http://chandra.harvard.edu/photo/chronological.htmli or X- 
ray binaries (microquasars) in our galaxy that have been ex- 
tensively studied (Guthmann et al.'2002' Fender'2004). On the 
other hand, the extremely relativistic particles that result from 
tidal acceleration outside the cone may interact with the ambi- 
ent medium within the Fermi system thereby transferring their 
tidal energy to the ambient particles that may escape from the 
system altogether and can appear as extremely energetic cos- 
mic rays. Such neutron star or microquasar particles may ac- 
count for the inferred ultrahigh energy primary cosmic rays 
with energies above lO^'^eV, since the corresponding extra- 
galactic particles would collide with the cosmic microwave 
background photons resulting in photopion production and pair 
creation (Greisen 1966 Zatsepin & Kuzmin 1966 1. It would 
therefore be interesting to search for any correlation between 
the directionality associated with the highest energy cosmic 
rays and the distribution of certain neutron stars and micro- 
quasars in our galaxy using the Pierre Auger Observatory. 
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